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We study the localization problem of one-dimensional interacting spinless fermions in an incom¬ 
mensurate optical lattice, which changes from an extended phase to a nonergoic many-body localized 
phase by increasing the strength of the incommensurate potential. We identify that there exists an 
intermediate regime before the system enters the many-body localized phase, in which both the 
localized and extended many-body states coexist, thus the system is divided into three different 
phases, which can be characterized by normalized participation ratios of the many-body eigenstates 
and distributions of natural orbitals of the corresponding one-particle density matrix. This is very 
different from its noninterating limit, in which all eigenstaes undergo a delocaliztion-localization 
transtion when the strength of the incommensurate potential exceeds a critical value. 

PACS numbers: 72.15.Rn,05.30.Fk,64.70.Tg 


Introduction.- The phenomenon of many-body local¬ 
ization (MBL) in interacting systems with disorder has 
attracted intensive studies recently Unlike the An¬ 
derson localization in noninteracting disordered systems, 
transition to MBL phase is not a thermodynamic phase 
transition in the general paradigm of phase transition. 
Instead, the MBL transition can be understood as a dy¬ 
namical transition from an ergodic phase to a nonergoic 
phase, where the thermodynamic equilibrium is not ac¬ 
cessible [:!. iTj, [Tot • While the MBL system differs from its 
noninteracting counterpart in dynamical features, such 
as dynamical correlations Q, the logarithmical growth of 
entanglement 11 T3, and conductivity 0,13 , the MBL 
transition can be also witnessed by analyzing the proper¬ 
ties of the system’s spectrum and many-body eigenstates 

B3- 

As most of previous works on disordered many- 
body systems focused on systems with random disor¬ 
ders, MBL in one-dimensional (ID) quasiperiodic sys¬ 
tems has also attracted great attentions recently [20l — 
3 | • This is not only due to the recent experimental ob¬ 
servation of MBL in interacting ultracold atomic gases 
trapped in quasiperiodic optical lattices [23j, but also 
the existence of some specific properties of quasiperiodic 
systems, e. g., t he presence of localization-delocalization 
transitions 
ID lattices 


24 , |25j and single particle mobility edges in 
261 ], which may bring new insights for under- 


standin g th e interplay of controllable disorders and inter¬ 
actions (21L [22J • In a recent work, Iyer et. al. presented 
numerical evidence for the existence of many-body er¬ 
godic and localized phases in an interacting qausiperiodic 
system [20| - however it is not clear whether an interme¬ 
diate phase exists between them, and if it exists, how to 
determine the boundary of different phases? Answering 
these questions is undoubtedly important for deepening 
our understanding MBL in the quasiperiodic system. 

Aiming to give answers to the above questions, in 


this work we study the interacting fermion model in a 
ID quasiperiodic optical lattice by using both the exact 
diagonalization method and density matrix renormaliza¬ 
tion group (DMRG) method. By analyzing the spec¬ 
trum and the properties of many-body wavefunctions, 
we identify that there exists an intermediate regime be¬ 
tween many-body ergodic and localized phases, in which 
the localized and extended states coexist. We can dis¬ 
tinguish the intermediate and MBL phases by analyzing 
the localization properties of the middle excited states 
and the energy level statistics. On the other hand, the 
boundary of extended and intermediate regimes can be 
fixed by studying the localization-delocalization transi¬ 
tion of the many-body ground state, which is found to 
be well characterized by the density distribution of the 
single-particle excitation of the interacting quasiperiodic 
system. 

Model and phase diagram.- We consider an interact¬ 
ing fermion model with the nearest-neighbor (NN) repul¬ 
sive interaction in a quasiperiodic optical lattice, which 
is described by 

H = Y^[~t(c]c j+ 1 + c] +1 £j) + hjhj + Vhjhj+i] ^ 

3 

with 

hj = hcos(2naj + 0), (2) 

where Cj is a fermionic annihilation operator, hj = CjCj 
the fermion number operator, h the strength of the 
incommensurate potential, a the irrational wavenum¬ 
ber, and 0 an arbitrary phase shift. The parameter V 
represents the strength of the NN interaction between 
fermions, which is permitted in a dipolar fermion sys¬ 
tem. For convenience, we set the hopping amplitude t to 
be the unit of the energy (t = 1), take a = v ^~ 1 , and 
consider the half-filling case, i.e., fixing N/L = 1/2 with 
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FIG. 1: Schematic phase diagram. In the Phase I and III, 
all of the many-body states are extended and localized, re¬ 
spectively. In the intermediate Phase II, the localized and 
extended states coexist. 
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FIG. 2: (Color online) The averaged 77 for (a) the ground state, 
(b) the mid one-sixth eigenstates, versus the quasidisorder 
strength h for half-filling systems with different sizes and fixed 
interaction strength V = 0.4. We use sample = 30, where a 
sample is specified by choosing an initial phase 6. 


N being the total number of fermions and L the lattice 
size. In the noninteracting limit with V = 0. the model 
reduces to the Aubry-Andre (AA) model [24(, for which 
all the single particle states are extended when h < 2 , 
but localized when h > 2. In the low-energy and weak- 
quasi-disorder limit, it was shown that the weakly in¬ 
teracting quasiperiodic system exhibits a metal-insulator 
transition by using bosonization technique and renormal¬ 
ization group analysis [271 ]. The MBL problem was only 
addressed very recently by studying the real-time dynam¬ 
ics of the interacting quasiperiodic system ( 20 | . 

In this work, we study carefully the transition pro¬ 
cess from the extended phase to the MBL phase. Be¬ 
fore going into the details, we present our main result 
firstly by displaying a schematic phase diagram in Fig. 
1 , which exhibits three phases, i.e., the delocalized phase 

I, the many-body localized phase III, and the intermedi¬ 
ate phase II. While all many-body states are extended or 
localized in the regime of I or III, respectively, localized 
and extended states coexist in the intermediate regime of 

II. In the following calculation, we shall fix V = 0.4 and 
determine the phase boundaries as marked in Fig.l. For 
the other values of V, the boundaries can be similarly 
determined. 

To characterize whether a given many-body state is 
localized, we consider the normalized participation ratio 
77 (NPR) m 


V(E) = 


1 


E{ ni ,n 2 ,...n L } hM{«i,n 2 ,...»iL})| V H ‘ 


( 3 ) 


which is a generalization of the participation ratio for a 
single-particle system, where i/)e ({%}) is the many-body 
wavefunction in the Fock basis with the eigenenergy E, 


and V[j is the Hilbert space dimension of the system with 
fixed total particle number ^2 = A’. In the thermody¬ 

namic limit, 77 tends to 0 for a focalized state, while it is 
finite for an extended state [22|. In Fig. [51 we display the 
change of 77 as a function of the quasi-disorder strength h 
by fixing the interaction strength V = 0.4 for the ground 
state and eigenstates in the middle regime of the energy 
spectrum under open boundary conditions. Although 77 
approaches 0 for both the ground state and middle ex¬ 
cited states in the large-h regime, there exists a regime 
around 2 < h < 2.4, in which 77 already approaches 0 for 
the ground state but remains finite for the middle excited 
states. In this regime, corresponding to the region of II 
in the proposed phase diagram, the many-body ground 
state is localized whereas the middle excited states are 
still extended states. 

To get an intuitive picture of localized and extended 
many-body states in the region of II, we next calculate 
the distribution of nature orbitals of the one-particle den¬ 
sity matrix 


p\4> a ) = n a \4> a ) (4) 

with the element of density matrix defined by 

p tj = {i>E\clcj\i/j E ), (5) 

where |t \>e) is a given many-body eigenstate, \<p a ) with 
a = 1,2,..., L are the natural orbitals, and the eigenval¬ 
ues n a are interpreted as occupations of natural orbitals 
which fulfills '^2 a n a = N. I 11 a recent work [lij], the 
one-particle density matrix has been applied to study 
the interaction-driven many-body localization in a dis¬ 
ordered fermion system and a step-like behavior in the 
occupation spectrum has been taken as an indicator of 
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FIG. 3: The distribution of n a for the ground state, the band 
center state and the highest excited state of the half-filled 
system with L = 16, V = 0.4, (a) h=2.1 and (b) h=3.0. We 
set ni > ri 2 > ... > ul, and use sample = 30. 
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many-body localization when we set n\ > 7r 2 > ■■■ > til. 
In Fig. (3[a) and (b), we display the distribution of n a 
(occupation spectrum) for the ground state, the highest 
excited state and the band center state of our model with 
the system size L = 16 in both the partially localized and 
many-body localization region, respectively. As shown in 
Fig. da) for the system with h = 2 . 1 , half of the natural 
orbitals for the ground state and highest excited state 
are almost fully occupied, i.e., (n a ) ~ 1 for a < N, 
and the other half are nearly unoccupied (n a ) ~ 0 for 
a > N, which indicates the ground state and highest ex¬ 
cited state are localized. However, for the center state 
in the band, the natural orbitals distribute continuously 
around the mean filling fraction n = N/L = 1/2 and 
no step-like change is detected, which suggests the band 
center state is still an extended state at h = 2.1. As a 
contrast, in Fig. [3])b), we show the occupation spectrum 
for the system at h = 3.0. It is obvious that a step-like 
jump in the occupation spectrum is developed for the 
band center state, which can be viewed as an indicator 
of many-body localization. 

To get a more comprehensive understanding, next 
we compare the energy resolved NPR.s in the three dif¬ 
ferent regions. In Fig. [4] we display the NPR versus the 
energy and its statistical histograms for systems with the 
same interaction strength V = 0.4 but three typical h , 
i.e., h = 1 , 2.2 and 2.8 from the top to bottom corre¬ 
sponding to the region of I, II and III in the proposed 


FIG. 4: Energy resolved NPR for the interacting fermion sys¬ 
tem with L = 18, N = 9 and V = 0.4. (a),(c) and (e) show 
r/(E) for all the eigenstates, (b),(d) and (f) show the corre¬ 
sponding histograms of NPR, for the system in the regime 
of I, II, and III, respectively. In (a) and (b), h = 1, all the 
states are extended; in (c) and (d), h = 2 . 2 , extended and 
localized states coexist; in (e) and (f), h = 2 . 8 , all the states 
are localized. 

phase diagram, respectively. From Fig. [2 (a) and (b), we 
see that 77 is finite in the whole permitted energy region 
when h = 1 , indicating all the states to be extended. 
On the other hand, when h = 2.8, r) approximately ap¬ 
proaches zero in the whole region as shown in Fig. [I] (e) 
and (f), suggesting that all the states are localized. How¬ 
ever, for the case of h = 2.2 in the middle region of II, 
the statistical histogram of rj exhibits quite different dis¬ 
tribution from the cases of h = 1 and h = 2 . 8 , i.e., both 
states with finite and nearly zero 77 are partially occu¬ 
pied, as shown in Fig. [2(c) and (d). In this region, the 
localized and extended many-body states coexist due to 
the interaction. To make a comparison, we also study 
NPRs of the corresponding non-interacting system in the 
same h-parameter regions and find no similar statistical 
histogram as that of Fig. [2 (d). 

To determine the phase boundary between the re¬ 
gion I and II more precisely, next we use DMRG method 






















4 


to study the ground state localization-delocalization 
transition. It is well known that the asymptotic behav¬ 
ior of a single-particle localized state is usually described 
by the exponential decay of its envelope. However, for 
a many-body localized state, no definition analogous to 
the exponential decay of a single-particle localized state 
can be constructed, as the many-body wave function de¬ 
pends on the coordinates of all N particles in the system. 
To characterize the many-body ground state localization- 
delocalization transition, instead we consider the density 
distribution of the single particle excitation, which is de¬ 
fined by 5rii = pn+i(i) — PN^i), he., the difference of den¬ 
sity distribution by adding one particle into the many- 
body system with Pn{i) = representing 

the density distribution of the ground state ip{N) with 
fixed N particles. Such a quantity has been successfully 
applied to characterize the localized edge excitation of 
many-body topological states 28, 29|. In Fig.[5ja) and 
(b), we display dn,; for the system with fixed V = 0.4 
and different h. For h = 1.9, the single particle excita¬ 
tion distributes over the whole lattice as shown in Fig. 
[5ja) , whereas it is localized in a narrow region for h = 2.1 
as shown in Fig.[5jb). This difference implies a signature 
of localization transition. To fix the transition point, we 

-32 


consider the inverse participation number [3C 
single-particle excitation, which is defined by 


IPR- 1 = 


of the 


( 6 ) 


with ]Th <5rii = 1. The IPR approaches 1 if the single¬ 
particle excitation is completely localized on a site. On 
the other hand, it tends to lattice size L if the single¬ 
particle excitation spreads over the whole lattice homoge¬ 
nously. In Fig.jS^c), we display IPR as a function of h, 
which decreases with the increase of h, indicating that the 
ground state wave function becomes more localized. For 
systems with different sizes, we observe that both IPRs 
have obvious drops around h = 2. We also display deriva¬ 
tives of the inverse participation numbers d(lPK)/dh in 
Fig.0d), which show obvious peaks in h = 2 , indicating 
the ground state localization transition point at h = 2 
for V = 0.4. 

Finally, we scrutinize the MBL of this system by 
using energy level statistics ii and the entanglement 
entropy 17( to determine the phase boundary between 
the region II and region III. A spectral distinction be¬ 
tween the many-body localized and the ergodic phases is 
based on the statistics of energy eigenvalues [3|, |4| . We 
use exact diagonalization to study gaps between adja¬ 
cent many-body levels under open boundary conditions, 
S n = E n+ 1 — E n ^ 0, where the eigenvalues of a given 
realization of the Hamiltonian E n are listed in ascend¬ 
ing order. Then we obtain the ratio of adjacent gaps as 
r n = min(S n , 5 n -i)/max(6 n , 5 n -i), and average this ra¬ 
tio over samples at each 0. For Poisson spectrum, the 
mean value of r is (r) ~ 0.387, a signature of a localized 



FIG. 5: (Color online) The distribution of single particle ex¬ 
citation for the half-filling system with L = 80, (a) h = 1.9 
and (b) h = 2.1. (c) IPR of the single particle excitation as a 
function of h. (d) The derivative of the inverse participation 
number d(IPR)/dh versus h. 


phase, and if the probability distribution of r is Gaussian 
orthogonal ensemble (GOE) random matrices, its mean 
value is (r) ss 0.529, a signature of a delocalized phase. 
As shown in Fig.^a) for the system with V = 0.4, the 
mean value (r) changes from 0.529 to 0.387, when h is 
increased from the extended region to the MBL region. 
We also study the entanglement entropy of the eigen¬ 
states averaged over the mid one-third states 17|. In 
Fig. ( 6 ](b), we display average entanglement entropy (S) 
and d(S)/dh as a function of h by fixing V = 0.4, where 
the entanglement entropy of a given eigenstate is defined 
as S = — JA A,7nA, ; , with Aj being the i-th eigenvalue 
of the corresponding reduced density matrix, which is 
obtained by tracing out one half of the system. While 
the average entanglement entropy follows a volume law 
in the extended phase, it decreases and eventually satu¬ 
rates at a constant independent of L deep in the localized 
phase with increasing the quasi-disorder strength, fulfill¬ 
ing an area law EmE Il5| . As shown in the inset of Fig. 
[ 6 ](b), the derivative of the average entanglement entropy 
presents an obvious divergent peak around the transition 
point, from which we estimate the MBL transition oc¬ 
curring at h c = 2.4 ~ 2.5, consistent with the result of 
NPR. 


Summary.- In summary, we have explored the many- 
body localization-delocalization transition of the ID 
quasiperiodic interacting fermion system and demon¬ 
strated that the phase diagram is composed of three dif¬ 
ferent phases, i.e., an extended phase, a MBL phase, and 
an intermediate phase. While all the many-body eigen¬ 
states are extended or localized in the extended or many- 
body localized phase, extended and localized states co¬ 
exist in the intermediate phase, which can be unveiled by 
the energy resolved distribution of NPRs and its statis- 
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FIG. 6: (Color online) (a) The ratio of adjacent energy gaps 
versus h. Here we use sample = 50 for L = 12 and L = 
14, sample = 30 for L = 16, and sample — 20 for L = 
18. (b) Averaged entanglement entropy ( S} and d(S)/dh as 
a function of h. Here we use sample = 500 for L = 8 and 
L = 10, sample = 100 for L = 12 and sample = 30 for 
L = 14. The interaction strength is fixed at V = 0.4. 

tical histograms. Furthermore, we have determined the 
phase boundary between the extended and intermediate 
phases by studying the density distribution of the single 
particle excitation and estimated the boundary between 
the intermediate and MBL phases by studying energy 
level statistics and the entanglement entropy. 
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